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Abstract 



We compare the properties of three-body systems obtained with two-body 
potentials with Pauli forbidden states and with the corresponding phase equiv- 
alent two-body potentials. In the first case the forbidden states are explicitly 
excluded in the calculation. Differences arise due to the off-shell properties 
of these on-shell equivalent potentials. We use the adiabatic hyperspherical 
method to formulate a practical prescription to exclude Pauli forbidden states 
• in three-body calculations. Schematic as well as realistic potentials are used. 

' Almost indistinguishable results are obtained. 



PACS number(s): 2I.60.Gx, 21.45.+V 



CI ■ 1 Introduction 



Cluster models are useful to describe spatially correlated systems when the intrinsic 
cluster structures essentially are frozen and the relative and intrinsic degrees of free- 
dom essentially are decoupled. A particular application during the last decade has 
been in studies of nuclear halos jjl |) and the concepts developed in this connection 
have already penetrated in to other subfields of physics || Q, ||. 

When different clusters contain identical particles the Pauli forbidden states 
must often be excluded. This is conceptually simple but often difficult in practise, 
since the proper antisymmetrization too easily spoils the advantages of the cluster 
description. For spatially extended systems with little overlap the details of ap- 
propriate approximations have so far been rather unimportant [^[ ^, |[ ^, [l0| Q |. 
However, beside the conceptual interest better treatments may be required by the 
need for higher accuracy. 

Furthermore, the validity of the cluster models may be extended to less pro- 
nounced halos by an accurate treatment of the exclusion principle. The dilemma 
is that the cluster models become more accurate with decreasing overlap while the 
overlap is necessary to provide the interesting structure and the binding energy. Too 
large an overlap destroys the fundamental decoupling of the intrinsic and relative 
cluster degrees of freedom, but the exclusion principle becomes important already 
at a smaller density overlap. 
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The correlations displayed in two and three-body halos have attracted a lot of 
interest and in particular those in light nuclei. If the halo consists of one nucleon 
bound to a nucleus the Pauli forbidden states are easily excluded by precise or- 
thogonalization. However, already three-body systems present more difficulties and 
approximations are often employed. Various prescriptions have so far been used in 
numerical calculations on halos consisting of two nucleons and a core || R . The 
simplest is to use a nucleon-core effective potential without bound states [H even 
when the nucleon then moves in a potential different from that of the identical nucle- 
ons within the core. The difference amounts in this way to an additional repulsion 
accounting for the Pauli exclusion principle. 

It is essential that this effective potential reproduces the nucleon-core low-energy 
scattering properties, i.e. scattering length and effective range or even better the 
phase shifts as function of energy. Therefore an obvious procedure jl2|, |l3| is to 
start with a deep potential appropriate for the nucleons in the core. Then construct 
a phase equivalent, approximate or precise, shallow or repulsive potential |14|j and 
compute the three-body observables of interest. The results could be compared to 
those obtained with the deep potential and the Pauli forbidden states excluded in 
one way or another [0, |l5) . 

This procedure must produce two sets of wave functions coinciding at large 
distances but differing at small distances by the additional node appearing inside 
the core by use of the deep potential. The idea is therefore that the short-distance 
behavior must be irrelevant when the cluster model provides a valid description. 
This conjecture was tested in three-body computations by using two potentials with 
the same scattering length and effective range [[zj and in two-body computations 
using completely phase equivalent potentials ||. 

The purpose of this paper is to investigate the consequences of using completely 
phase equivalent two-body potentials || [l^, [l3[ [hi], [ll| in three-body calculations 
based on the adiabatic hyperspherical approach, see e.g. [[| ||. The method was 
suggested and previously tested approximately for two neutrons and a core j?J . In 
this report we shall carefully follow this prescription using phase equivalent poten- 
tials and rigorously derive the short and long-distance behavior of the all-decisive 
effective radial potentials. We shall consider cases when one, two or three of the 
two-body potentials have Pauli forbidden states characterized by various quantum 
numbers. 

We begin with a brief sketch of the general method in section 2. We then in 
section 3 give the characteristic properties of the constructed phase equivalent two- 
body potentials, where we use a simple s-wave potential as an analytical example. 
In section 4 we derive the short and long-distance properties of the three-body 
adiabatic effective radial potentials. In section 5 we show a number of qualita- 
tively different numerical examples. Finally section 6 contains a summary and the 
conclusions. 

2 Adiabatic hyperspherical expansion 

In the adiabatic hyperspherical method the coordinates of the three-body system 
are divided into the slow (hyperradius) and fast (angular) variables. The eigenvalues 
and eigenfunctions for the fast angular degrees of freedom are then calculated for 
fixed values of the slow radial variable. These eigenfunctions arc then used as a basis 
for expansion of the total wave function. Eigenvalues for the fast angular variables 
serve as adiabatic potentials for the slow radial variable. The great advantage of 
this basis is that it correctly reproduces the long-range asymptotic behavior of the 
three-body systems including the pathological Efimov effect. 
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2.1 Hyperspherical coordinates 



We first introduce the usual Jacobi coordinates |J for three particles with masses 
rrij and coordinates r,: 



1 rrijmk . 1 mi(rrij + rrik) ( jnjYj+m^r^ 



y m TOj + mfc V m M \ mj + nik ) 

where M = mi + mi + and m is an arbitrary mass scale. The hyperspherical 
coordinates now consist of a hyperradius p = yxf + yf and five dimensionless 
angles Q — {c^, fl X i, Q y i}, where f2 X j and £l V i specify each of the directions of 
and yi and on = arctan(xi/?/i). The kinetic energy operator T can be written in 
terms of the hyperspherical coordinates as 

T = ^+2^' < 2 > 

9 2m P dp 2p + 4 2mp 2 ' [ ' 

Id 2 f 2 T 2 

A n = - • /o n 5 ? sin ( 2a ) ^ 4 + T-f- + ~ , ( 4 ) 
" sm(2a) 9a 2 v ' sin 2 a cos 2 a w 

where Z| and ^ are the angular momentum operators related to the x and y coor- 
dinates. 

2.2 Adiabatic expansion 

The three-body hamiltonian is given by 

H = T p + ^Al+Y / V t ( P ,n), (5) 

where Vi is the interaction between the two particles j and fc. Here {i,j, k} is 
a permutation of {1,2,3}. The eigenvalue problem is solved in two steps. First 
the slow coordinate p is fixed and the (always discrete) spectrum X n (p) and the 
corresponding angular eigenfunctions (p, fi) are calculated as function of p. For 
this one can either use the Schrodinger equation 

Un + £ v * ~ A «(-°)) ^ °) = > ( 6 ) 
or the equivalent Faddeev equations 

(^- A " ( ^tR2^y + ^- V;<i, " ( ' 9 ^ ) = °'^W> ( 7 ) 

(i) 

where the Schrodinger wave function <f> n is the sum of the Faddeev components </>„ 
or more precisely $ n (p, H) = £V <^(p, 0)/ sin(2aj). 

The total wave function ^(p, f2) is expanded on the complete angular basis set 
i.e. 

n ni 

where the phase space factors p~ 5 / 2 and [sin(2ai)]~ for convenience are extracted 
explicitly. 
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The radial expansion coefficients f n (p) are then solutions to 

'' "' { UP) + - Qnn(p) + v 3 ( P ) E ] f n {p) 



2m dp 2 2mp 2 \ 4 
d_ 
d~p 



2Pnn>0>)lT-Qnn>(p))fn>(p), H') 



= E 

n' T^n 

where V3 (p) is an anticipated three-body potential for fine tuning and the coupling 
terms P and Q are given as 

p m Ap) = J <m$* n (p,n)—$ nl (p,n), (10) 

Qnn'W = J cM<i>* n (P,n)-^<l> n ,(p,n) (11) 

with the angular volume element defined by dfi = sin 2 a cos 2 adadtt x dQ y . When 
all but the lowest terms are neglected we obtain the usual adiabatic approximation. 

The angular eigenvalues A„ (p) are therefore closely related to the adiabatic radial 
potentials. The behavior of these angular eigenvalues decisively determines the 
properties of the three-body system. We shall therefore analytically and numerically 
investigate X n {p) for phase equivalent two-body potentials differing by the number 
of bound states. We shall especially study the asymptotic properties analytically. 



3 Phase equivalent two-body potentials 

Let us assume that a two-body system with relative coordinate r, reduced mass p 
and potential V has at least one bound state with the energy E\ S j — —h 2 K 2 /2/u, and 
the total wave function ipi S j(r) — ui S j(r)r~ 1 yi S j m expressed in terms of the radial 
wave function ui S j(r) and the spin-angular part yi S j m , where (s,l,j) are the spin, 
the orbital and total angular momenta, respectively, while m here is the projection 
of j on the z-axis. We then construct a new potential 

V pep (r) = V(r) + AV lsj (r) , (12) 
Wsj(r) = -PjPiPs j-^lnH \u lsj (r')\ 2 dr'^j , (13) 

where (Pj, Pi, P s ) are projection operators on states with the corresponding angular 
momentum quantum numbers. Thus the additional potential AVi S j (r) is an effective 
radial potential which only is active for the spin and angular momentum quantum 
numbers of the bound state in question. The potential V pep (r) is only missing 
the bound state ipi S j of V(r), but otherwise these two potentials have identical 
properties, i.e. they are phase equivalent with the same bound state spectrum 

If desired this construction can be repeated to eliminate possible other bound 
states of the new potential. The series of potentials then only differ by the number 
of bound states. 

The wave function ui S j (r) has the following asymptotic forms at short and large 
distances 

uisj(r -> 0) cx r l+1 , u Uj (r ->■ 00) oc e~ Kr . (14) 

The additional term AVi S j in Eq. (fL2|) d iverges as r~ 2 at small distances and vanishes 
exponentially at large distances ]9|, |12|, |l3| , |i~4| , |l5[ ] 

fi 2 91 3 

AV lsj (r^O)^ PjPtPs P , (15) 

P 7* 
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4frV 

AV lsj (r -► oo) -> PjPiPs exp(-2 K r) . (16) 

A* 

We shall in the remaining part of the present paper often consider an original 
potential V with one bound s-state. This is the most important case in the practical 
applications for halo states and a generalization is easy and straightforward by use 
of Eqs.(|l^) and (13). Properties and procedure can be studied in more details in 



simple cases. For illustration we therefore construct an analytical s-wave potential 
and the corresponding solution labeled by "0". The wave function Uo(r) must be 
proportional to r at small distances and must fall off exponentially at large distances. 
One of the simplest functions obeying these constraints is 

Uo(r) — sinh(/3r) exp(— ar) , (17) 

which is the solution with the energy Eq for the potential Vb given by 

Eo = ~(a - (3) 2 = , V (r) = ^-2/3a [1 - coth(/3r)] , (18) 

where K = a — (3 > 0. 

The asymptotic behavior of this potential is then 

4-2 / 2cv \ 2?i 2 

V (r -> 0) = — + 2(3a + 0(r) , V (r -» oo) = /3aexp(-2/3r) , (19) 

2/i V t J M 

where 0(r) is a function vanishing at least as fast as r when r decreases towards 
zero. 

The phase equivalent potential in Eq.(|l2|) has the additional term AVq defined 
in Eq.(|l3|). For the example in Eqs.jT^) and ( p^|) we then get 

AVo(r) = — 8a(3 2 (a 2 - (3 2 ) sinh(/3r) 
M 

a(e 2ar + 1) sinh(/3r) - f3(e 2ar - 1) cosh(/3r) 
(a 2 (cosh(2/3r) - 1) - (3 2 (e 2ar - 1) + a/3 sinh(2/3r)) 2 ~ 
h 2 16a/3 2 (a 2 — (3 2 )e ar sinh(/3r) [acosh(ar) sinh(/3r) — /3sinh(ar) cosh(/3r)] 
M [a 2 (cosh(2^r) - 1) - ^ 2 (e 2Qr - 1) + a/3 sinh(2/3r)] 2 

which only acts on s-waves and asymptotically behaves in agreement with Eqs.(|l5| 
and (PI 



4 Asymptotics of the angular eigenvalues 

The different two-body potentials from section 3 should now be used in the three- 
body computations described in section 2. Two-body interactions are needed for 
the three two-body subsystems. Each of these can in principle support Pauli for- 
bidden states of specific but different properties. In any case the angular eigenvalue 
spectrum, or equivalently the adiabatic radial potentials, is crucial for the behavior 
of the resulting three-body system. We shall therefore calculate, as far as possible 
analytically, the leading terms in the asymptotic expansion of the adiabatic ra- 
dial potentials both for small and large p. Numerical examples for all p-values are 
postponed to the next section. 
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4.1 Short distances 



Let us assume that the original two-body potentials are of short range and finite 
at r = 0. The different angular momentum states then decouple in the three-body 
Schrodinger equation in the small-distance limit p = 0. For each set of quantum 
numbers l x and l y referring to the i'th set of Jacobi coordinates, the angular eigen- 
value equation in Eq.(|^) reduces to 

1 52 .sin(2a 4 ) ^ + ■ ^ + 1 > 



sin(2aj) da 2 1 sin 2 a. 



2mp 2 



^^•UnW H =^n Ws > (21) 



where A„ = v 2 — 4 and the potentials are kept for later use although their contri- 
butions vanish for p = 0. The eigenvalues and the corresponding eigenfunctions are 

1,0, Pi 

v n = 2(n + l) , n = 0,1,2,... (22) 
®mj y = ^ko P n* +imv+1/2) (cos 2on)Y lmmm (n xi )Y lymy (Slyi) , (23) 

where Yi m and p^ +1 ^ 2 ' l+v1 ^ 2 ' are t ne spherical harmonics and the associated Leg- 
endre polynomials, respectively. 

Let us now assume that only one bound state of angular momentum l x is Pauli 
forbidden in one of the two-body subsystems, which we choose to relate to the x- 
coordinate and denote by i. The phase equivalent potential in Eq.(|l2|) then has an 
additional term AV^, which for small p contributes to Eq.(^l"[) as 

2m P\ V p ~ 2m P 2 ( 3 + 2 ^) fi2 S 2(3 + 2Uj3 

— — t\V fi m -j - fi - — fx , (24) 

h h pr z sm oti 

where we used that pr 2 /m = p 2 sin 2 a%. All other partial waves remain unchanged. 

The Schrodinger equation in Eq.(|2l]) is then for this particular l x modified at 
small p into 

O.m n 2 . — . \ ~ 

(25) 

where we omitted the index i on the coordinates and tildes indicate the new eigenval- 
ues and eigenfunctions. Thus the Schrodinger equation is changed as if the angular 
momentum quantum number l x were increased by two units. For all other partial 
waves with l x different from the quantum number of the two-body bound state the 
Schrodinger equation remains unchanged. The new spectrum and the corresponding 
eigenfunctions are given by 

V n = 2{n + 2) , n = 0,1,2,... (26) 

= iR5S) Pt +5,2 ' ly+m (cos 2a)Y lxmx (Q x )Y lymy (Q y ) . (27) 

Thus, except for the missing lowest eigenvalue, the three-body angular spectrum 
for the phase equivalent potential is identical to the original potential up to the first 
two terms in the expansion for small p, i.e. 

, 2 = (2(n+l)) 2 + ^]>>(0) , £ 2 = (2 (n + 2)) 2 + ^5>(Q), (28) 
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where n = 0, 1, 2, .... 

The next step is to include more than one two-body bound state in the same 
subsystem. If the orbital angular momenta of these states are different the three- 
body angular equations decouple and each of them is analogous or identical to 
Eq.(|2^). Thus each Pauli forbidden state is removed independently. On the other 
hand if the angular momenta are equal the procedure for one two-body bound 
state is repeated and two units have to be added to the corresponding l x both in 
Eq.(|2^) and in the Legendre polynomial of Eq. (|2lj) , i.e. p^ +b / 2 ' l y +1 / 2 ^ should be 
substituted by p( l *+ 9 / 2 ' l y+ 1 / 2 *) wn ile the spherical harmonics remain unchanged. 
The spectrum is again identical to that of the original potential except for the 
missing two lowest eigenvalues. 

Let us now assume that more than one of the two-body subsystems have Pauli 
forbidden states. Using the phase equivalent potentials the three-body Schrodinger 
equation for small p then becomes 



1 d 2 . ,„ N T 2 2 , ^2(3 + 2$) 



sin(2a0 + - Jf- + — §— + V V , 2 xj ' P. 



(0) 



sin(2a l ) da 2 % sin 2 a, cos 2 a, ^-f sin 2 ay h 

$„l = V^nL , (29) 




where j in the first summation includes subsystems with a Pauli forbidden bound 
state of angular momentum l^j . The projection operator projects on states with 

the angular momentum l x = in the j'th Jacobi coordinate system. When more 
than one subsystem has a Pauli forbidden state the different partial waves couple 
already in the limit p = 0. Therefore we only maintained the conserved total angular 
momentum L as an index on the solutions. 

The spectrum arising from Eq.(|2|) can not be obtained analytically. Let us then 
first repeate that treating each of the additional potentials AV while neglecting the 
others results in Eq. (^7|) . Thus simultaneous inclusion of more than one of these po- 
tentials related to different subsystems could be expected to move the spectrum by 
approximately the sum of that obtained by each of the individual treatments. The 
corresponding approximate wave functions are analogously obtained from Eq.(p7j) 
with the appropriate addition of units to the angular momentum quantum numbers 
in the Legendre polynomials. However, this is not the exact solution due to the cou- 
pling of different partial waves. If the states with l x °^ and are Pauli forbidden 
in the subsystems labeled "1" and "2", respectively, all l x \ in subsystem "1" couple 
provided they receive finite contributions from rotation of the state with l x °2 ■ Thus 
numerical solutions are required to get the precise solutions. 

A better, still approximate, analytically obtained spectrum can be found by 
first solving the equation when l x \ differs from all three l x °^ and all additional 
AV terms are neglected. The wave functions are then the usual hyperspherical 
harmonics. Then when l x \ equals one and only one of the values l^j we use the 
wave functions in Eq. (p7|) . When l x \ is equal to two or three identical -values 

-,(i x +9/2,Z H + l/2) 



the Legendre polynomial in Eq.(|27D has to be substituted by P,i' xTJ/ *- Ly+±/ r 



p^(2 x +i3/2,2„+i/2) ^ reS p ec ^j ve |y The spectrum is finally found approximately by using 
these approximate wave functions $^ ; and first order perturbation AE n i x i y on top 
of the hyperspherical spectrum where one or more of the lowest levels are missing. 
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The energy correction is then 




£<*■ 



3 




^niO > ( 30 ) 



where again the summation is over the Pauli forbidden states. This means that at 
most two non-zero terms are present. Unfortunately it not easier to use this approx- 
imation than to compute the spectrum numerically without any approximations. 

In conclusion, in the limit of small p the original and the phase equivalent poten- 
tials have, except for the Pauli forbidden eigenvalues, exactly the same spectrum of 
angular eigenvalues when all Pauli forbidden states occur in one two-body subsys- 
tem. When more than one subsystem support Pauli forbidden states the spectrum 
is only approximately the same for p = 0, again excluding the lowest-lying Pauli 
forbidden levels. 

4.2 Large distances 

The hyperspherical spectrum is for short-range two-body potentials also obtained in 
the limit of p = oo as established previously j|, (Tt|, [l8| . It is essentially equivalent 
to the spectrum at p — given by Eq.(p2f). However, in addition some of the angular 
eigenvalues diverge parabolically towards — oo in a one-to-one correspondence with 
the bound states in the two-body subsystems. The corresponding wave functions 
are the hyperspherical harmonics in Eq. (p3|) except for the lowest diverging levels 
where the structure at p = oo is that of two particles in a bound state while the 
third particle is in a continuum state relative to this two-body structure Q . 

The convergence towards the hyperspherical spectrum is faster the higher the 
angular momentum quantum numbers. The limits are reached at distances com- 
parable to the effective ranges of the two-body interactions [|[ [l0| with the crucial 
exception of s-waves, where the two-body scattering lengths define the distance for 
convergence |]lq| . The different partial waves decouple in the sense that contribu- 
tions to a given state from higher angular momenta are at least one order smaller 
in expansions in terms of p^ 1 . Thus for example p- waves only receive contributions 
of the same order in the p _1 expansion from the s-waves, which in turn only receive 
contributions of the same order from other s-waves. These couplings are basically 
determined by the centrifugal barriers increasing with angular momentum. 

The relatively fast convergence for the higher angular momentum states isolate 
the s-waves as the only interesting contributors at large distances. The diverging 
levels are of special interest, since some of them correspond to Pauli forbidden states, 
which must disappear by changing to the phase equivalent potentials in Eq. . 
In fact the goal of the present subsection is to compare the three-body angular 
eigenvalues of deep and phase equivalent two-body potentials differing only by the 
number of bound states. This should then allow us to design a prescription for 
removing Pauli forbidden states at large distance. 

The necessary subsequent connection from small to large distances requires 
knowledge of the approach to the asymptotic behavior, i.e. at least a minimum 
of information about rate of convergence and wave function decomposition during 
the approach. Here the s-waves are crucial and we shall therefore derive a num- 
ber of characteristic features at large distances, i.e. corresponding to distances p 
simultaneously outside the effective ranges of all the two-body interactions. These 
properties are revealed by use of the zero-range potentials with the related analytical 
solutions. 

For short-range potentials and large distances the Faddeev equations are more 
convenient than the Schrodinger equation. We therefore start with Eq.(|7]) for s- 
waves, which means that the partial angular momenta l x = l y = in all three Jacobi 
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coordinates. We choose the system labeled i and transform the other two Faddeev 
components into this system, project out the s-wave components and multiply by 
sin(2ai). The wave function resulting from <!>„ after these operations is denoted by 
$n and from Eq.(|^) we then obtain the equation of motion as 



d 2 
da 2 



^\a i ) + ^-V i (J^p S ina i )^(a i )=0, i = 1,2,3, (31) 



$W(a)=«A«(a)+]T 



in 
1 



sm(2ipij) J w 



where tfij 



da'^ia'), (32) 
and P{ijk) is the 



(— l) p ( 4 -7 fe ) arctan ^-\/ m fc( m i + m 2 + m 3 ) / (mirrij) 
parity of the permutation {ijk} while m and pt = mjmk/(mj + mk) are the nor- 
malization mass m and the reduced mass, respectively. The last term in Eq.(|32"|) 
arises from the transformation of the two Faddeev components different from i, see 
for example || ^|, [ll| . A complete solution for square well potentials can be found 
ing§. " 

For large distances and short-range potentials Vj vanish when on is much larger 
than the largest effective range of the potentials divided by p. In this region the free 
solution is obtained for each Faddeev component q>n . In the remaining vanishingly 
small region of c^-space around = the potential Vi contributes. In the region 
of small on the integral in the function in Eq. ( p2| ) can be expanded in terms of 
on. To first order we then obtain 



$«(«<<!) = </>«(«) + 2a ^ 



sin(2</jy ) 



i = 1,2,3 . 



(33) 



For large distances p ^> R = max(i?i, i?2, R3), where Ri is the range of the potential 
Vi, the solution to Eq.(|3l|) in the a^-region, scj = y / ^psina i > R, is the free 
solution 

sin v (a - - ) , (34) 



where Aj is constant. The structure of this solution is by definition independent 

of the potentials. The wave function <f>„ (Vij) m Eq. (|33|) is then given by Eq. (p4|) , 
since the argument ip^ is finite. 

Increasing the parameter p in Eq.(|3l|) effectively decreases the range and in- 
creases the strength of the potential and thus the zero-range limit is approached. 
The related eigenvalue problem can then be reformulated in the limit p — > 00 for 
zero-range potentials. For x$ < Ri and on -C 1 Eq.([3l|) is approximately given as 



27:^ 2^7) r, ( ,-,)<i»:,"( ( ,) = o ^ /= 1.2.:; 



h 2 - 2 



(35) 



and for Xi > Ri the potential vanishes. The Schrodinger equation for then 
arises by adding the three equations in Eq.(^5|). 

The singular points are on = or Xi = and the boundary conditions for these 
zero-range potentials are then expressed in terms of the two-body scattering lengths 

di as 



1 d 



L<I> 



W dx, 



1 



1 d 



W da. 



(i) 



a i= 



1,2,3, 



(36) 



where wc define the modified scattering lengths by 2j = dj ^/ pi/m. Thus in the zero- 
range limit the eigenvalue spectrum is determined entirely by the scattering lengths 



9 



and identical spectra results from potentials with the same scattering lengths. In 
particular the original deep potential and its phase equivalent partner must then 
produce identical angular spectra at large distances, where the zero-range limit is 
reached. 

The boundary conditions in Eq. d36| ) with the explicit expressions for the angular 
functions in Eqs.(^) and (Q) result then in a system of linear equations for the 
coefficients A, i.e. 



AiV cos I 



This system of equations has non-trivial solutions only when the corresponding 
determinant D is zero, i.e. D = det{My} = with 



M = 2 Bin[,(^-f)] 
sm(2cpij) 



j 7T\ p . / 7T\ 

i>cos v— + — sin v— 
V 2 J 5i \ 2/ 



. (38) 



When two or three of the scattering lengths are infinitely large D = has a specific 
negative solution v which results in the so-called Efimov effect [jlTj . 

For finite scattering lengths there are two types of solutions to D = in the 
large-distance limit where p » max(|oi|). One type of solution is only present when 
there is a bound state with binding energy Ej, in one of the two-body subsystems 
j. This is most easily seen if we assume that aj < and \dj\ 3> Rj which imply 
that the binding energy Eb « h 2 /{2maf) = h 2 /{2maf}. Then D = has a solution 
v which diverges linearly with p, i.e. 

P \ 2 a ~ 2m P 2 TP ^ 'qo^ 
v = —i— , A = v -4si o— i/b • (39) 

Thus this angular eigenvalue diverges parabolically with p at large distance. 

The other type of solutions has eigenvalues v approaching constants as p — > oo. 
This is understood from the dominance of the terms proportional to p in Eq. ( |38| ) 
resulting in the solutions 



Const 



, sin (iS) = , v n = 2(n + 1) , n = 0, 1, ... . (40) 



displaying a spectrum identical to that obtained for small p. These eigenvalues v 
can be found to next order in an expansion in powers of p by expanding D as 

£> = det{M^} = p 3 ^^sin 3 (VJ) (41) 



' «^ ^os(,-)sm 



where 0(/o) at most increases as p for p — > oo. These leading two terms then give 
the eigenvalue equation D = 0, i.e. 



sin 



("~) + v cos -J2 Oi = , (42) 



P 

which has an explicit solution for the inverse function p(v), i.e. 
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The eigenvalue solutions v n to Eq. (|42|) are then to first order in p 1 given by 




(44) 



7T p 



A„ = ^-4= ( 2 ( n+l))*-4-^±^^, 

7T p 



(45) 



where the leading terms already are given in Eq. (|40|) . 

In conclusion, the phase equivalent potential has at large distances the eigen- 
value spectrum in Eqs.(^) and ( (iB]) while the original deep potential has the same 
spectrum plus the additional eigenvalue in Eq.(|39|) resulting from the bound state. 
More two-body bound states would simply result in precisely the same number of 
parabolically diverging eigenvalues of the form in Eq.(|4(]). Therefore except for 
this (or these) additional eigenvalue(s) the angular spectra of the two potentials are 
asymptotically identical in the large-distance limit. 

4.3 Connecting small and large-distence eigenvalues 

The original and the phase equivalent potentials have identical angular eigenvalue 
spectra at both large and small distances when only one two-body subsystem has 
bound states. This agreement includes at least the first two terms in the expansion 
both for small and large p. The deep potential has in a one-to-one correspondence 
with the two-body bound states in addition a number of low-lying eigenvalues, which 
diverge parabolically towards — oo for p — > oo. The structure of the related states is 
asymptotically precisely that of the bound two-body states with the third particle 
in the continuum. When more than one two-body subsystem has bound states the 
small-distance behavior of the two spectra is only approximately identical. 

We can now formulate a convincing prescription for dealing with Pauli forbidden 
states in three-body cluster computations. In complete analogy with the two-body 
problem, where the lowest eigenvalues of the two-body hamiltonian are forbidden, 
we declare the lowest angular eigenvalues, or equivalently the lowest adiabatic po- 
tentials, as Pauli forbidden. We consequently remove them from the basis for all 
values of p in the subsequent solutions of the radial equations. The resulting spectra 
for the original potential with the Pauli forbidden states and its phase equivalent 
partner are then approximately identical for all values of the hyperradius. 

The small differences must be compensated by slightly different three-body in- 
teractions, see Eq.(^). This adjustment is necessary in any case to ensure that the 
three-body binding energy is sufficently precise. Even two-body interactions repro- 
ducing the two-body phase shifts would otherwise underestimate the three-body 
binding energy. The three-body potential is related to polarization of the parti- 
cles beyond that included already in the two-body potentials. It must then arise 
when all three particles interact simultaneously, i.e. are close to each other. The 
three-body potential may also be viewed as an adjustment of the off-shell behavior 
of the two-body potentials exclusively determined by the on-shell behavior. After 
the fine tuning practically indistinguishable and almost exactly identical results are 
obtained by use of the two different two-body potentials with different appropriate 
three-body potentials. 

5 Numerical illustration 

The three-body angular eigenvalue spectra are identical or at least approximately 
identical for short distances for phase equivalent two-body potentials. In addition, 
the large-distance behavior is determined by the scattering lengths and therefore 
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Figure 1: The lowest angular eigenvalues A n as functions of hyperradius for the 
original schematic potential and its phase equivalent partner, see Eqs.(|l8|) and (po|). 
The parameters are a = 4.151ml- 1 and = 1.063 £m _1 , h 2 /{2(i) = 41.4 MeV fm 2 
and 77i 3 = 18/i. The other two potentials are equal and given by V13 = V23 = —7.8 
MeV exp(— (r/2.55fm) 2 ). The lowest parabolically diverging level corresponding 
to the Pauli forbidden state has been removed. Only s-waves are included in the 
computation. 



also identical for such potentials. The intermediate distances must be studied nu- 
merically. 

We shall first consider a three-body system where only one two-body subsystem 



has a bound state. We choose the schematic potentials in Eqs.(18) and (|20|) for 
these two particles. The parameters are chosen to reproduce the neutron-neutron 
scattering length of 18.45 fm and the effective range of 2.83 fm. The two remaining 
interactions with the third particle are assumed to be equal and given by V13 — 
V23 = —7.8 MeV exp(— (r/2.55fm) 2 ) ||. The masses correspond to two neutrons 
and 9 Li. The three-body system then resembles the halo nucleus n Li ( 9 Li+n+n) 
except that now the deep neutron-neutron potential has a bound state. 

The resulting lowest angular eigenvalues are shown in Fig. [l] when only s-waves 
are included in the calculation. The lowest diverging eigenvalue of the original 
deep potential with the Pauli forbidden state is removed. The two-body potentials 
behave as r _1 at the origin and the eigenvalues are therefore according to Eq.(p8|) 
linear functions of p for small distances around zero. The hyperspherical spectrum 
is obtained at large distances. The two spectra are almost indistinguishable for all 
distances from zero to infinity. 

When more than one two-body subsystem has bound states the short-distance 
behavior of the angular spectra for phase equivalent potentials are only approxi- 
mately equal. We show in Fig. || numerical examples for realistic s-wave potentials 
for three-body calculations of 9 Li+n+n, where we assume spin-zero for the 9 Li- 
core. With three-body interactions for fine tuning these potentials all produce a 
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Figure 2: The lowest three-body angular eigenvalues A n as functions of hyperra- 
dius for realistic s-wave potentials for the 9 Li+n+n system. The neutron-neutron 
interaction is from 

0- We 

use three different neutron- 9 Li potential with the same 
low-energy scattering properties, i.e. (i) a deep potential with Pauli forbidden 
states, V nc — —288 MeV exp(— (r/1.15fm) ), (ii) its phase equivalent partner with- 
out bound states and (hi) a shallow potential without bound states, V nc = —5.6 
MeV exp(— (r/2.55fm) 2 ) from pj] |. Only levels antisymmetric with respect to in- 
terchange of the two neutrons are shown. We also removed the lowest parabolically 
diverging level corresponding to the Pauli forbidden state. Only s-waves are in- 
cluded in the computation. 



11 Li structure in agreement with measurements. By construction the s-wave phase 
shifts for the deep and phase equivalent potentials are identical for all energies, and 
indistinguishible from those of the shallow potential for energies below 15 MeV. 

The two spectra arising from the truly phase equivalent potentials, (i) and (ii) 
in Fig. |^, are again almost indistinguishable although they differ slightly more at 
small and intermediate distances than those in Fig. |l|. Otherwise the features of 
the schematic potential remain for these two potentials. On the other hand the 
spectrum from potential (iii) differs qualitatively from the other two spectra in 
spite of the fact that the low-lying s-wave phase shifts are identical for all these 
three potentials. The hyperspherical spectrum is still for the shallow potential 
obtained in both limits of p = and p — oo, but now the lowest level starts at 
A = like the (removed) Pauli forbidden level for the deep potential. The phase 
equivalent potential does not have any level at A = 0. The large-distance behavior 
coincide precisely for all three potentials. The angular eigenvalues are in general 
more smooth for the shallow potential. 

The different behavior at short distance must of course be reflected in different 
phase shifts at higher energies, but low-energy properties are only marginally in- 
fluenced. For the same reason the spatially extended halo structure resulting from 
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Figure 3: The same as Fig. H but now we included both s and p-waves in 
the neutron- 9 Li interaction. The Pauli forbidden states are two-body s-waves 
and the interactions for s-waves are as in Fig. 2. For p- waves we use the 
interaction in ]TT[ denoted potential IV, with a central interaction given by 
—5.0 MeV exp(— (r/2.55fm) 2 ), and a spin-orbit interaction given by 33.6 MeV 
l„ c • s„ exp(— (r/2.55fm) 2 ). These interactions are realistic for use in calculations of 
halo properties. 



these three potentials can hardly be distinguished. This may also be expressed as 
a dominance of intermediate and large distances in the radial wave function. 

The s-wave potentials used in Fig. || are realistic choices for computations of 
the three-body properties of 11 Li. The s-waves are dominating for this system, but 
the contributions from p-waves are also necessary. These partial waves couple in 
the three-body computation and several new angular eigenvalues appear as seen in 
Fig. |^. Comparing to Fig. ^ we can identify the new levels which predominantly 
are of p-wave character. Even with these couplings the quantitative features of the 
three potentials are still remarkably similar with the exception of the lowest level 
for the repulsive s-wave potential. 

We also want to compare the spectra when the p-waves dominate in the three- 
body structure while the Pauli forbidden state still remains in one of the two-body 
s-wave potentials. A halo nucleus with these features is 6 He ( 4 Hc+n+n) where 
the shallow s-wave potential for n Li is substituted by a realistic repulsion. All 
three potentials (deep, phase equivalent, repulsive) are adjusted to reproduce the 
measured low-energy s-wave phase shifts. All three angular eigenvalue spectra are 
shown in Fig. ^ when only s-waves are included. 

The repulsion is revealed at small distance by the increasing levels which rel- 
atively fast join the corresponding levels from the other two potentials. For the 
shallow s-wave potential in Fig. |2j the levels from the two phase equivalent po- 
tentials varied in contrast relatively fast at small distances, but the outcome is in 
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Figure 4: The same as Fig. g for the 4 He+n+n system. The deep potential is given 
by —75.06 MeV exp(— (r/1.43fm) 2 ) and now the shallow potential is substituted 
by the repulsive potential 48 MeV exp(— (r/2.33fm) 2 ). The measured low-energy 
s-wave phase shifts are reproduced by all potentials. 



both cases that all spectra quickly becomes very similar and at large distances com- 
pletely indistinguishable. The overall quantitative agreement between these spectra 
is again remarkably good. 

Including s, p and d- waves in the computation produce more levels and a number 
of avoided crossings as seen in Fig. [|. However, the larger variation with hyperradius 
does not spoil the close quantitative resemblance of the three spectra. Instead the 
agreement is actually improved due to the coupling of the different partial waves 
and the related avoided crossings. For example the s-wave increasing from the 
starting point of A = 32 on Fig. [| is now on Fig. || quickly joining another level also 
originating from 32 but with higher angular momentum. Thus, the angular spectra 
are also remarkably similar when the p-waves dominate as in realistic computations 
of the three-body properties of the 6 He system. 

The radial wave function is obtained from Eq.([)|) and the lowest is often carrying 
more than 95 % of the probability. To demonstrate the similarity of the results found 
by use of the different potentials it is therefore sufficient to show the lowest function 
fi, see Fig. ^. Indeed the quantitative agreement is overwhelming and the curves 
are essentially indistinguishable except perhaps at the very small hyperradii (blown 
up and shown in the inset). The slightly smaller probability around the hyperradius 
of 1 fm for the repulsive potential simply reflects that the effective potential is larger 
in this region, see Fig. ||. On the other hand the probability is already small and 
the difference is therefore insignificant. Even these qualitatively different two-body 
potentials lead to very similar radial wave functions provided the phase shifts are 
identical. The two phase equivalent potentials produce practically the same radial 
wave functions for all hyperradii. 
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Figure 5: The same as Fig. | but now we included both s, p and <i-waves. The 
s-potentials are as in Fig. ||, while the p, d and spin-orbit interactions are taken 
from [ fro[ . The measured low-energy s-wave phase shifts for all partial waves are 
reproduced by all potentials. The Pauli forbidden states are two-body s-waves. 

6 Summary and conclusion 

We use the adiabatic expansion and the hyperspherical coordinates to solve the 
three-body problem for given short-range two-body potentials. When the three par- 
ticles are clusters consisting of identical fermions the antisymmetrization between 
clusters is technically difficult and a number of different practical approximations 
are used. For nuclear halos the particles are often nucleons interacting with a nu- 
cleus. Then this problem can be formulated as exclusion of Pauli forbidden states 
within the nucleus from the space suspending the physical three-body states. In 
two-body calculations this is straightforward by direct projection but the use of 
phase equivalent potentials still provide an interesting alternative. 

We use phase equivalent two-body potentials with a different number of bound 
states and compare the corresponding three-body angular eigenvalue spectra. We 
show analytically that these spectra, or equivalently the adiabatic radial potentials, 
at large distances are identical to second order, except for a number of diverg- 
ing eigenvalues precisely corresponding to the additional Pauli forbidden two-body 
bound states. This limit is mathematically obtained by use of zero-range potentials 
and the complete analytical solution is therefore derived and discussed. 

We show furthermore analytically that the spectra at small distances, with pre- 
cisely the same exception for the additional bound states, also are identical to second 
order provided only one of the two-body potentials differ in the number of bound 
states. When more than one of the two-body potentials differ in the number of 
bound states the spectra at small distances are only approximately identical. Wc 
give an estimate in perturbation theory of this relatively small deviation. 
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Figure 6: The lowest radial wave function as function of hyperradius for the 
4 He+n+n system for the potentials used in Fig. ||. The inset shows the small 
distances where the wave functions differ from zero and that of the repulsive poten- 
tial from the two other potentials. 

In the numerical illustrations we use schematic as well as realistic potentials 
for the two-neutron halo nuclei n Li ( 9 Li+n+n) and 6 He ( 4 He+n+n). The bound 
states are for convenience assumed to be s-states and the two halo nuclei then 
provide tests on systems where respectively the s and p-configurations dominate 
in the three-body structure. The three-body angular spectra for phase equivalent 
potentials are in all cases in remarkable numerical agreement also at intermediate 
distances. The radial wave function is finally computed for the different potentials. 
They turn out to be almost indistinguishable for strictly phase equivalent potentials. 

However, we find a small difference at very small distances when we use an 
approximate two-body potential reproducing only the low-energy phase shifts with- 
out complete phase equivalence. Such phcnomcnological potentials have often been 
used in practical computations of nuclear halo properties. Due to the large spatial 
extension of halos the probability is by definition very small at small distances and 
the deviation is therefore insignificant for all practical purposes. 

We formulate a prescription to deal with the Pauli principle for two-nucleon 
halos. Assume that the two-body potentials have Pauli forbidden states. We then 
construct and use the phase equivalent partners without these forbidden states. 
The difference between the three-body adiabatic potentials arising from the origi- 
nal and the phase equivalent two-body potentials is essentially only additional basis 
functions in the hyperspherical adiabatic expansion. Removing these from the adi- 
abatic basis makes the resulting adiabatic spectra asymptotically identical both at 
small and large distances and approximately equal as well at all the intermediate 
distances. 

The removed basis functions correspond at large distances precisely to two par- 
ticles in Pauli forbidden states and the third particle in a continuum state. At 
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small distances the lowest levels of the original deep potentials correspond to iden- 
tical fermions occupying the same states. Removal of these terms therefore forces 
the particles to occupy higher-lying orbits and thereby introducing the necessary 
repulsion preventing violation of the Pauli principle. In conclusion, this method to 
exclude the Pauli forbidden states in a three-body system has firm mathematical 
and numerical foundations. It is a practical and accurate alternative to the other 
existing methods. 
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